Two Transitions in the Damping of a Unitary Fermi Gas 
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We measure the temperature dependence of the radial breathing mode in an optically trapped, 
strongly-interacting Fermi gas of 6 Li, just above the center of a broad Feshbach resonance. The fre- 
quency remains close to the unitary hydrodynamic value, while the damping rate reveals transitions 
at two well-separated temperatures, consistent with the existence of atom pairs above a superfluid 
transition. 



Optically-trapped, unitary Fermi gases 0, 0] test pre- 
dictions for exotic systems, from nuclear matter [j. Q 
and quark-gluon plasmas 6] to high temperature super- 
conductors 7]. Recent measurements of the heat capac- 
ity of a unitary gas reveal a transition 0, . This has 
been interpreted as the onset of superfluidity, using a 
Bardeen-Cooper-Schrieffer (BCS) - Bose-Einstein Con- 
densate (BEC) crossover approach which was initially de- 
veloped for high-temperature superconductors jsL 1 1 0I| . In 
a unitary Fermi gas, pair interactions between particles 
are "strong" in the sense that the zero-energy scatter- 
ing length is much greater than the interparticle spac- 
ing, as achieved by tuning near a Feshbach resonance . 
Such a gas exhibits universal features Q, 0, 0] . At suf- 
ficiently low temperatures, unitary Fermi gases are be- 
lieved to comprise normal atoms, condensed pairs and 
noncondensed pairs @, 0|. Fermionic atom pairs are 
probed in recent projection experiments Il2l Il3l Il 1| and 
in measurements of the pairing gap fl5l Il6|. Evidence 
for superfluid hydrodynamics in a unitary gas appears 
in anisotropic expansion and in the breathing mode 
frequencies and damping rates 0, 0, 0] . 

In this Letter, we report the precision measurement of 
the temperature dependence of the frequency and damp- 
ing rate for the radial breathing mode of a unitary Fermi 
gas of 6 Li. We identify two transitions in the damping 
rate which occur at well-separated temperatures. This 
is consistent with a qualitative picture of a unitary gas 
in which the superfluid transition temperature lies well 
below the temperature at which noncondensed pairs first 
form 0, Ho| . Remarkably, neither transition is accompa- 
nied by an abrupt change of the frequency, which in the 
whole range of temperatures remains within a few per- 
cent of the unitary hydrodynamic value. Below the lower 
transition temperature, the damping rate extrapolates to 
zero at zero temperature, as expected for a superfluid. 
The higher temperature transition can be interpreted as 
arising from the breaking of noncondensed atom pairs by 
the collective excitation. 

In the experiments, we prepare a degenerate 50-50 mix- 
ture of the two lowest spin states of 6 Li atoms b y fo rced 
evaporation Q in an ultrastable CO2 laser trap |2fJ ] . At 
a bias magnetic field B of 840 G, just above the center 
of the Feshbach resonance [U , the trap depth is low- 
ered by a factor of ~ 580 in a few seconds Ui 17] and 



then recompressed to 4.6% (for most of the experiments) 
of the full trap depth in 1.0 s and held for 0.5 s to assure 
equilibrium. A controlled amount of energy is added to 
the gas by releasing the atoms from the trap for a short 
time and then recapturing the cloud 0, 0| . The gas is 
then allowed to thermalize for 0.1 s. 

The radial breathing mode is excited by releasing the 
cloud and recapturing the atoms after 25 /Lis (for 4.6% 
trap depth). After the excitation, we let the cloud oscil- 
late for a variable time thoid, at the end of which thegas is 
released and imaged after ~ 1 ms of expansion 0, • 

Radial breathing mode frequencies uj and damping 
times r are determined from the oscillatory dependence 
of the released cloud size on thoid 0, E El- For 
each temperature, 60-90 values of thoid are chosen in the 
time range of interest. These values of thoid are ran- 
domly ordered during data acquisition to avoid system- 
atic error. Three full sequences are obtained and aver- 
aged. The averaged data is fit with a damped sinusoid 
xq + Acxp(— t/r) sin(cjt + ip). We have obtained oscilla- 
tion curves at 30 different temperatures, containing data 
from 6300 repetitions of the experiment. 

For most of the data reported, the total number of 
atoms is N = 2.0(0.2) x 10 5 . From the measured trap 
frequencies, corrected for anharmonicity, we obtain for 
4.6% trap depth: ujj_ = y/uj x uj y = 2tt x 1696(10) Hz, 
Lu x /uj y = 1.107(0.004), and uj z = 2ir x 71(3) Hz, so 
that ui — (ojxLUyUJz) 1 ^ = 2tt x 589(5) Hz is the mean 
oscillation frequency and A — lu z /uj± = 0.045 is the 
anisotropy parameter. The typical Fermi temperature 
Tp = (3NY^ 3 hu>/kB of a corresponding noninteracting 
gas is ~ 2.4 /iK, small compared to the final trap depth 
of U /k B = 35 yK (at 4.6% of full depth). The coupling 
parameter of the strongly-interacting gas at B = 840 G 
is fcpa ~ —30.0, where fikp = \J2m ksTp is the Fermi 
momentum, and a — a(B) is the zero-energy scattering 
length estimated from Ref. j2l| . 

The dimensionless empirical temperature T is deter- 
mined by the method implemented in 0, : The column 
density of the cloud is spatially integrated in the axial di- 
rection to yield a normalized (integrates to 1), one dimen- 
sional, transverse spatial distribution n[x). This distri- 
bution is fit to determine the empirical reduced temper- 
ature T using a finite temperature Thomas-Fermi profile 
with a fixed Fermi radius, which is measured in a sepa- 
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rate experiment at the lowest temperatures [8J, |9J ■ The 
empirical temperature T is numerically calibrated to the 
theoretical reduced temperature T/Tp 0,0]. In Ref. 0, 
we show that a simple approximation relating T to T/Tp 
is given by, 



The first method for estimating (p 4 )/{p 2 ) assumes 
that the spatial distribution of the gas is nearly a zero- 
temperature Thomas-Fermi profile, which is a good ap- 
proximation at the lowest temperatures. In this case, the 
corrected frequency is [ijj Er| 



t fV tt? 



(1) 



Eq. [2 yields accurate values of T/T F for f > 0.45 
and provides a reasonable estimate at lower temper- 
atures, where higher precision can be obtained using 
the calibration |9J[. Here (3 is the unitary gas parame- 
ter P, 0, IE EE Hj , which we recently measured to be 
fj = -0.49(0.04) (statistical error only) 0,01. 

The mode frequency provides important information 
on the state of the system. The frequency versus the em- 
pirical temperature for experiments at 4.6% trap depth is 
shown in Fig.^ The figure shows the measured frequen- 
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FIG. 1: Frequency u) versus empirical reduced tempera- 
ture T: Raw data (empty circles with error bars); Data 
corrected for anharmonicity: Black dots- using a zero- 
temperature Thomas-Fermi profile and Triangles- using a 
finite-temperature Thomas-Fermi profile. The dot-dashed 
line is the unitary hydrodynamic frequency ujh = W10/3u>±. 
The dashed line at the top of the scale is the frequency 2to x 
observed for a noninteracting gas at the lowest temperatures. 

cies Wmeas (open circles), uncorrected for anharmonicity 
in the trapping potential, as well as the frequencies after 
correction by two different methods. The fre quen cy cor- 
rection is proportional to the ratio (p 4 ) / (p 2 ) |25j, where 
p is the transverse radius of the expanded cloud. For an 
isentropic unitary gas, we obtain |26j 
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with bn a hydrodynamic scale factor [lj, |27| ■ For 1 ms of 
expansion, we obtain bn — 13.3. 
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where (x 2 ) is the transverse mean square width of the 
gas in the x-direction after the expansion. Applying this 
method over the whole temperature range, we find that 
the corrected frequencies (black dots in Fig. ^) remain 
very close to the unitary, hydrodynamic value, shown as 
a dot-dashed line. 

In the second method, we include the effects of the 
finite temperature on the spatial profile of the cloud 
by calculating the ratio (p 4 )/(p 2 ) = (4/3) (x 4 )/(x 2 ) di- 
rectly from the fitted one-dimensional finite temperature 
Thomas-Fermi profiles. The corrected frequencies are 
displayed in Fig. ^ as open triangles. The higher or- 
der corrections have been calculated and are found to be 
negligible. 

The radial breathing frequency varies smoothly over 
the whole temperature range, and remains close to the 
value loh = yl0/3cjj_ = 1.83 u) ± predicted by hydro- 
dynamic theory for a unitary gas, where l/(kpa) = 
[H HI US iH HI H3. Such temperature indepen- 
dence has been observed previously in a BEC F° r 
the unitary Fermi gas, the observed frequencies are far 
from 2uj x — 2.10 oj±, the value observed for a noninter- 
acting gas at the lowest temperatures, which is shown as 
a dashed line at the top of Fig.^ This, as well as the ob- 
served hydrodynamic expansion, justifies the use of the 
hydrodynamic expansion factor bu in Eqs. Eland over 
the whole temperature range. 

The frequencies obtained using the finite temperature 
corrections (triangles in Fig. ^) rise 4% above u>h at the 
highest temperatures. This can be explained by the slow 
decrease in the collision rate of the unitary gas at higher 
temperatures j^, which makes the gas slightly less hy- 
drodynamic and pulls the frequency up toward the non- 
interacting gas value. The slow increase in frequency 
is consistent with our previous estimate of the reduced 
temperature for ballistic expansion of the unitary gas, 
T/T F > 3 

In contrast to the frequency, the damping rate of the 
radial breathing mode reveals two transitions, at T ~ 0.5 
and at T ~ 1.0, as shown in Fig. El 

The lower transition at T — 0.5 is indicative of a super- 
fluid phase transition, consistent with our previous study 
of damping versus temperature |l7j . Between f = 0.1 
and T — 0.5, the damping rate decreases with decreas- 
ing T, the opposite of the behavior expected for a de- 
generate, collisionally hydrodynamic gas pjrj l3?l l3S| and 
consistent with a superfluid picture [Tj| • This behavior is 
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FIG. 2: Temperature dependence of the damping rate for 
the radial breathing mode of a trapped 6 Li gas at 840 G, 
showing two transitions. Solid dots are the main data set 
taken at 4.6% trap depth. Other symbols are for the system 
with scaled parameters: Two squares - at 0.85% of full trap 
depth; Four diamonds - at 19% of full trap depth; Two open 
circles - at 3 times smaller number of atoms. The dot-dashed 
line is the maximum damping rate for a classical harmonically 
trapped gas with binary collisions. The solid line is Eq. @ 
which extrapolates close to zero at zero temperature. 

similar to that observed in the radial breathing mode of a 
BEC 0. For 0.65 <f< 1.0, the damping rate is nearly 
independent of temperature. This change in behavior is 
accompanied by a significant, reproducible notch in the 
damping rate near T = 0.6. Using Eq. ^ we find that 
T = 0.5 corresponds to T/Tp = 0.35. This transition 
temperature is somewhat higher than the predicted su- 
perfluid transition temperature of T c = 0.29 Tp, as well 
as the value T c = 0.27 Tp estimated from the observed 
slope change in the heat capacity, after temperature cal- 
ibration |9j. However, it is not clear that the observed 
change in the damping rate should occur precisely at the 
same temperature as for the change in heat capacity 0- 

Damping below T = 0.5 may arise from the interaction 
between a superfluid core and a normal gas in the edges 
of the cloud. Normal fermionic excitations are present 
even when the core is superfluid [3{|, because the local 
Fermi energy at the edges can be smaller than /cbT c . It is 
possible that the transition near T = 0.5 arises as the core 
of the cloud changes from superfluid to collisional as the 
temperature increases. Indeed, in the temperature region 
from 0.65 < T < 1.0, the damping rate is close to the 
maximum allowed for a collisional gas 1/r = 0.09 uj± |37| . 

The higher temperature transition region, 1.0 < T < 
1.2, may arise from the breaking of noncondensed pairs. 
In this region, the frequency remains near the uni- 
tary hydrodynamic value while the damping rate rises 
roughly linearly with increasing T from approximately 
the maximum value allowed for a collisional gas, 1/r = 
0.09 oj± [37j, to a value roughly 1.5 times larger, signal- 



ing the appearance of a new channel of energy loss. The 
temperature range over which this transition occurs is 
much smaller than the temperature scale, of order Tp, 
over which the binary collision rate of a unitary gas de- 
creases in this temperature region j3f|. Although the 
frequency is nearly constant, this increase in damping 
rate resembles that which we observed at low tempera- 
ture, at a magnetic field of 1080 G, where the coupling is 
reduced to kpa = —1.35 [19j. Similar behavior was first 
observed in Ref. 0] and attributed to a pair breaking 
process The enhanced damping in all of these ex- 

periments may be of the same nature: both the decrease 
in coupling and the increase in temperature reduce the 
pairing gap, making it comparable to the collective exci- 
tation energy Tluj and therefore, causing fermion pairs to 
break. 

The second transition region occurs at high temper- 
atures where the energy of a strongly-interacting Fermi 
gas m erges with that of an ideal noninteracting Fermi 
gas |8j, |9j, |l0| . Using Eq.^ we see that the second transi- 
tion region 1.0 < f < 1.2, corresponds to 0.71 < T /Tp < 
0.86, close to the temperature range estimated for the 
vanishing of noncondensed pairs [9j, |lfj, |40| . In principle, 
merging of the ideal and unitary gas energies near Tp can 
arise (at least in part) because the unitary gas becomes 
classical at high temperatures. However, the appearance 
of enhanced damping supports the concept that the bind- 
ing energy of noncondensed pairs is decreasing as part of 
the merging process |9j, llOj , causing collective excitations 
to break pairs. 

We have investigated the possibility that the observed 
variation of the damping rate with temperature might 
arise in part from oscillations of different components 
of the gas at different frequencies. In a Bose-Einstein 
condensate with a thermal cloud, such behavior leads to 
revivals of the net oscillation amplitude, altering the ap- 
parent decay rates |4lj . In the present experiments, we 
find no evidence for such revivals, even after increasing 
the time over which the decay of the mode is observed. 

The near linear dependence of the damping rate on T 
in the region 0.1 <T< 0.5 is well fit by 

— = 0.146 (0.004) T - 0.0015 (0.0014), (4) 

r uji_ 

for the main data set taken at 4.6% trap depth and N — 
2 x 10 5 atoms. The damping rate extrapolates close to 
zero at zero temperature, consistent with a superfluid. 

We have examined the dependence of the damping rate 
on the trap oscillation frequency lu± and on the number 
of atoms N. Dimensional analysis requires that 1/r = 
uj± f(T/Tp,N, A), where / is a dimensionless function. 

For fixed T/Tp (or fixed T), we find that the function / 
cannot have a strong number dependence. For example, 
it cannot be oc UbTf / {%uj±_) cx iV 1 / 3 or its inverse. This 
is established by examining the scaling of the damping 
rate with the atom number. We find that the damping 
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rate at 4.6% of full trap depth does not change apprecia- 
bly when the number is reduced by a factor of ~ 3: The 
damping rate at reduced number, open circles in Fig. |3 
lies very close to the main data set (solid dots) when 
plotted versus T. Hence, it is likely that 1/r depends on 
N only via the combination T/Tp, and the most general 
formula for 1 jr is then limited to 

-=^/(^,A). (5) 

T Tp 

Experimentally, we are not able to test whether the 
damping rate depends on A. We have verified that 
1/r versus T scales approximately as by monitor- 
ing the breathing mode in the trap at 0.85% of full 
depth (uj± = 728(4) Hz , squares in Fig. 0) and at 19% 
(u>± = 3343(20) Hz, diamonds). In both cases, 1/r in 
units of lo±_ is comparable to that of the main data set. 

Quantum viscosity ^3 recently has been suggested 
as the mechanism for the small damping rate observed 
for the axial breathing mode in Ref. |l8(. The quan- 
tum viscosity r\ is of order hkp/cr, where the collision 
cross section a oc l/k F in the unitary limit []J. Hence, 
•q oc hkp oc Tin, where n is the density. For our system, 
the radial damping rate arising from quantum viscosity 
is estimated to be 1/(twj_) = 3 x 10~ 5 . This is consis- 
tent with the extrapolated value at T = 0, but cannot 
explain the observed rates of order 0.014 w± at our low- 
est temperature. Hence, the decay of the radial mode is 
probably by a different mechanism. 
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